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Abstract 



In the framework of spherical geometry for jelHum and local spin den- 
sity approximation, we have obtained the equilibrium values, fs{N,Q, of 
neutral and singly ionized "generic" A^-electron clusters for their various spin 
polarizations, C- Our results reveal that fs{N,(^) as a function of C behaves 
differently depending on whether corresponds to a closed-shell or an open- 
shell cluster. That is, for a closed-shell one, fs{N, Q is an increasing function 
of C over the whole range < ^ ^ 1) and for an open-shell one, it has a 
decreasing part corresponding to the range < C ^ Co) where Co is a po- 
larization that the cluster assumes in a configuration consistent with Hund's 
first rule. In the context of the stabilized spin-polarized jellium model, our 
calculations based on these equilibrium values, fs{N, (), show that instead 
of the maximum spin compensation (MSG) rule, Hund's first rule governs the 
minimum-energy configuration. We therefore conclude that the increasing 
behavior of the equilibrium values over the whole range of is a necessary 
condition for obtaining the MSG rule for the minimum-energy configuration; 
and the only way to end up with an increasing behavior over the whole range 
of C is to break the spherical geometry of the jellium background. This is the 
reason why the results based on simple jellium with spheroidal or ellipsoidal 
geometries show up MSG rule. 



I. INTRODUCTION 

Since the production and study of sodium clusters by Knight et a/.,0 the physics of metal 
clusters has attracted much interest .iH Metal clusters are composed of atoms and have 
properties that are different from both a single atom and the bulk metal. However, when 
increasing the size of the cluster, its properties evolve to those of the bulk. The many-body 
technique suitable for these systems is the density functional theory (DFT). It is a well- 
known fact that the properties of alkali metals are dominantly determined by the delocalized 
valence electrons. In these metals, the pseudopotentials of the ions do not significantly affect 
the electronic structure, because the Fermi wavelengths of the valence electrons become much 
larger than the metal lattice constants if one replaces the ions with a uniform positive charge 
background. This fact allows us to replace the discrete ionic structure by a homogeneous 
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positive charge background. This approximation is known as the jelhum model (JM). The 
simplest way of applying the JM to metal clusters is to replace the ions of an A^-atom 
cluster by a sphere of uniform positive charge density and radius R = (zNY^^rg, where z 
is the valence of the atom and is the bulk value of the Wigner-Seitz radius of the metal. 
However, since the ionic density near the surface of a metal differs from that of the bulk 
region, one may resort to the diffuse jellium model (dif-JM)i in which the density of the 
jellium background falls to zero, from the bulk value, within a length of a few atomic sizes. 
Application of the dif-JM to metal clusters results in a better agreement of the theory and 
experiment over the JM results. However, in spite of its simplicity and success in predicting 
some properties of bulk metals and metal clusters, the JM, which was originally developed 
for bulk metals, has some drawbacks.i0 In order to overcome the deficiencies of the JM, one 
should take some details of the ionic structure into account. Among the various methods of 
improvement, the first attempts that kept the simplicity of the JM and overcame some of the 
deficiencies of the JM resulted in the development of the stabilizedjellium model (SJM)0 
or pseudojellium model.llilli The SJM was applied to metal clustersBll and improved some 
results of the simple JM. Montag et al.^ in their structure-averaged jellium model (SAJM), 
which added the ionic surface energy to the SJM energy functional, made it more suitable 
for metal clusters. 

In other approaches, some researchers relax the spherical geometry and use the JM with 
spheroidal or ellipsoidal shapes,EHli which are suitable for open-shell clusters. Relaxing 
the shape of the jellium background as well as its density distribution while keeping charge 
neutrality at evervDoint in space, called the ultimate jellium model (UJM), was introduced 
by Koskinen et alB^ Since the bulk density in the UJM (r^ = 4.18) is close to that of sodium, 
its results can be compared with experiments on Na metal clusters. However, since the 
jellium density and the electron density in the UJM are locally equal everywhere, the UJM 
cannot describe the ionized clusters. 

In a recent work, by taking the spin degrees of freedom into account in the process of 
stabilization, we have generalized the SJM to the stabilized spin-polarized jellium model 
(SSPJM).i In the SJM, the equilibrium bulk density is a free parameter and the experi- 
mental value is used for it. However, in the SSPJM the equilibrium bulk density parameter, 
fs{oo, C), is polarization dependent and to the best of our knowledge, no experimental data 
are available for it. Therefore in the SSPJM we take 



C(oo, C) = rt(oo, 0) + Ar^(oo, C), (1) 

in which f^(oo, 0) is the equilibrium bulk value for the spin-compensated system which takes 
the experimental value of metal X, and Ar^*^(oo, () is obtained by the application of the 
local spin-density approximation (LSDA) to the infinite electron-gas system. All equations 
throughout this paper are expressed in Rydberg atomic units. It turns out that Ar^'^(oo, () 
is an increasing function of (. By taking ( = {N^ — Ni)/N for a monovalent metal cluster, 
and calculating the corresponding Ar^^(cxo,C) , it is possible to determine the appropriate 
radius for the spherical jellium through R{N,() = N^^^fs{oo,C). Therefore, any variation 
in N-^ and iV|, keeping N = N-^ + constant, leads to different values for ( and thereby, 
different values for -R(A^, C)- 

Application of this model for the energy calculation of metal clusters show that the energy 
is minimized for a configuration with maximum spin-compensation (MSC).ii That is, for 
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clusters with an even number of electrons, A^, the number of up-spin electrons, A^|, equals 
the number of down-spin electrons, N^; and for an odd number of electrons, — = 1. 
This MSG rule leads to the fine structure in A2(A^) (see Fig. 5 of Ref. ^3|) and the odd-even 



alternations in the ionization energies [Fig. 7(c) of Ref. |23 



In the present paper, keeping the spherical geometry for the jellium, we have inves- 
tigated the finite-size effects on the equilibrium values, and their consequences on the 
SSPJM calculations. First, using the LSDA, we have found fs{N,(), the equilibrium 
values of the closed-shell A^-electron neutral and singly ionized "generic" clusters with 

= 2,8, 18,20,34,40 for all possible polarizations. By generic we mean no specific metal 
but a simple jellium sphere which can assume its equilibrium size for any given value of A^. 
By fitting the calculated results to a polynomial with even powers of ( (in the absence of 
magnetic fields, the physical properties are invariant under the transformation ( 
we have found two analytic equations for Ars{(), one for neutral and the other for singly 
ionized clusters. Employing these analytic equations in the SSPJM calculations for Na 
clusters shows that as before, the MSG rule is at work and the results do not show any 
significant changes over our previous results on the ionization energies.ii In the next step, 
we have found the values fs{N, () for all different neutral and singly ionized generic clusters 
(A^ < 42) with different spin configurations {0 < ( < 1), and have shown that fs{N,() 
behaves differently for open-shell and closed-shell clusters. That is, we have found that for 
a closed-shell cluster (and its two nearest neighbors) it is an increasing function of ( over 
the whole interval < C ^ 1? whereas for an open-shell cluster (except for the two nearest 
neighbors of a closed-shell cluster), it has a decreasing behavior over < ( < Co, and an 
increasing behavior over Co < C ^ 1- However, in both open-shell and closed-shell clusters, 
the global minimum of energy, i.e., the ground-state energy, corresponds to a configura- 
tion in which fs{N,() is a minimum. Here, (q corresponds to an electronic configuration 
for which Hund's first rule is satisfied. By subtracting the ground-state energy of a singly 
ionized generic cluster from that of its neutral counterpart, we obtain the ionization energy 
of that generic cluster. Galculation of these ionization energies for A^ < 42 shows a good 
agreement with experimental results on Na clusters. We see that, in the ionization-energy 
plot, although the saw-toothed behavior remains, the pronounced shell effects near closed 
shells, seen in the simple JM results, are substantially reduced. Thanks to the appreciable 
reduction of the shell effects in the ionization energy results of the above-mentioned calcu- 
lations, we have performed the SSPJM calculations using the set of values fs{N,()- The 
results show that, instead of the MSG rule, Hund's first rule is governing the ground-state 
configuration. We have also performed simple JM-LSDA calculations for the ground-state 
energies of Na clusters but, instead of using the ordinary bulk value (3.99), we have used 
the increasing function fs{oo,C) as in Eq- (0)- The results show that, here also, Hund's 
first rule remains at work and therefore, we conclude that with spherical geometries, the 
increasing behavior of Arf*^(cxD,C) is a necessary condition (but not sufficient) for realizing 
the MSG rule, and the energy corrections over the simple JM due to stabilization are also 
needed. Therefore, in order to improve the SSPJM results one should insist on the increasing 
behavior for fs{N,()- The only way which guarantees such behavior for all values of A^ is 
to relax the spherical constraint on the geometry and consider ellipsoidal shapes for open 
shells. Using the LSDA and ellipsoidal geometry, one could obtain fs{NX), for a given set 
of values of A^ and C, by finding the values of the ellipsoid axes which correspond to the 
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minimum of energy. We expect that in this case fs{N, () will be an increasing function of ( 
for both open-shell and closed-shell clusters. Then, if one performs the SSPJM calculations 
for ellipsoidal clusters using these new increasing functions fs{N,(), one would obtain the 
MSG rule (and thereby the odd-even alternations) with improved results. 

The organization of this paper is as follows. Section y is devoted to calculational schemes. 



In Sec. Ill we present the results of our calculations. In Sec. IV we conclude the work 



II. CALCULATIONAL SCHEMES 

The energy functional in the SSPJM is given byil 

EsspjM[n^,ni,n+\ = EjmK, n|,n+] + (^mI^) +WR{n)) J dr n+(r) 

+ {SvUs{n) J dr 0(r)[n(r) - n+(r)], (2) 

in which Ejm is the energy fuctional of the simple JM, is the Madelung energy, wr is the 
average value of the repulsive part of the pseudopotential, and (5f )ws is the average of the 
difference potential over the Wigner-Seitz cell and the difference potential, 6v, is defined as 
the difference between the pseudopotential of a lattice of ions and the electrostatic potential 
of the jellium background. n|(r) and nj(r) are, respectively, the up-spin and down-spin 
electron densities with the total electron density given by n{r) = n|(r) -|- ni{r); and n is 
the uniform jellium density which is ^-dependent; for our previous SSPJM calculations we 
have used the values given by Eq. (|^). In the spherical JM, the functions will depend only 
on the radial variable. The function B(r) becomes a simple radial step function, 6{R — r), 
where R = N^^^fs with = (S/Anny^^. In the case__pf jellium with a sharp boundary, 
n+(r) = n6{R — r), whereas for the diffuse case we tal- 



, t n{l-{R + t)e--f^/*[sinh(r/t)]/r}, r <R , . 

~ 1 n{l - {{R + t)/2R){l - e-2^/*)}i?e(^-^)/7r, r > R, 



where t is a parameter related to the surface thickness. 

To evaluate the total energy of a cluster, we solve the Kohn-Sham (KS) equations self- 
consistently. The effective potential in the KS equations for the SSPJM calculations is given 
by 



where 



<ff(K.'^i.^+];r) = r) + <^([n-f,n|];r) + ((5i;)ws W0(r), (4) 



Kn,l;r) = 2/<(r'M£l^^. (5) 



which appears in the electrostatic part of the total energy of the simple JM energy functional 

+ 2/ dr (f){[n,n+];r)[n{v) ~ n+{v)], (6) 
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and 

ClK^^-i];!-) = (y=]d- (7) 

For E^f. we use the LSDA with the Perdew-Wang parametrization for the correlation 
part.B 



III. RESULTS AND DISCUSSIONS 

In the first step, by solving the KS equations for spherical geometries of the jellium in the 
LSDA, and finding the minimum energies, we have obtained fs(iV, the equilibrium val- 
ues of the closed-shell neutral and singly ionized jellium clusters with = 2, 8, 18, 20, 34, 40 
electrons. In the calculations we have considered all different possible polarizations. That 
is, for an A^-electron cluster (A^ = N^ + Ni), we have found fs{N, Q for all different polariza- 
tions (0 < C < 1) corresponding to the configurations Ni = N^,Ni = N^ — l,Ni = N^— 2, 
. . . , A^l = 0. The results are shown in Figs. |l|(a) and 0(b). We see that for these closed-shell 
clusters, fg is an increasing function of C,. By a least-square fitting of the results to the 
polynomial 

TsiO = ao + aaC' + a^C + "eC", (8) 

we have obtained = 4.28, 02 = 2.15, 04 = —2.41, og = 1.84 for neutral, and = 4.62, 
02 = 1.49, 04 = —3.31, ag = 3.92 for singly ionized clusters. In Figs. |l](a) and |l](b), we have 
also shown plots using Eq. for the two sets of the coefficients and compared with the 
bulk function f^*^(cx), Q. The values of show that the simple JM with LSDA predicts a 
larger atomic spacing for a cluster than the bulk value (oq > 4.18). Regardless of the Oq 
values (since we are studying the variations of with respect to C), we have performed our 
SSPJM calculations using 

Ar,(C)=a2C' + a4C' + a6C', (9) 

with corresponding coefficients for neutral and singly ionized clusters. Figure |^ compares 
the plots of Eq. (^ for neutral and singly ionized clusters with Arf^(cxD,C) of the bulk 
jellium. We see that they are increasing functions of C,. Using Eq. (^), we have performed 
the dif-SSPJM calculations for Na clusters with fs(oo,0) = 3.99. After the self-consistent 
calculations of the KS equations, we have obtained the dif-SSPJM ground-state total energies 
of the neutral and singly ionized Na clusters (A^ < 42). The calculations show that the MSG 
rule is governing the ground-state configuration. In the dif-SSPJM calculations, we have 
taken t = 1. In Fig. ^ we have plotted the ionization energies and compared them with 
our previous resultscJ and also with the experimental values. As is shown, there are no 
significant changes over the previous results. 

In the second step, using JM-LSDA we have found the values fs{N,Q for all different 
neutral and singly ionized jellium clusters (A^ < 42) with different spin polarizations. Here, 
for a cluster with specified values of A^ and (, the value fs(A^, () minimizes the total energy 
of the cluster. In Fig. |^ we have shown the values of fs{N,() and corresponding total 
energies per electron, E{N, C)/N, as functions of the number of the electrons, A^, for neutral 
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jellium clusters. The two plots show the same structure. That is, the maxima and minima of 
these two plots correspond to the same values of A^. The values of the ground-state energies 
correspond to polarizations Co consistent with Hund's first rule and form the Hund curve. 
In each of the plots, the uppermost value for a given N corresponds to the configuration 
with MSG. The values in between correspond to intermediate polarizations. We therefore 
conclude that for an open-shell cluster the functions fs{N,() and EJN.()/N have minima 
for a polarization Co 7^ which is consistent with Hund's first rule.E3Ea That is, for open- 
shell clusters (except for the nearest neighbors to the closed-shell cluster) and E{N, () are 
decreasing functions of ( for < ( < (q and increasing functions for Co ^ C ^ 1- the bulk 
jellium, the function f^'~'(oo,C) is an increasing function over the whole range < C ^ 1- 
In Fig. H we have compared the functions fs(27, C) and fs(34, C) with ff'^(oo,C)- We see 
that for the closed-shell cluster (A^ = 34) the volume is an increasing function of C on the 
whole range < C ^ 1- But for the open-shell cluster (A^ = 27) the volume decreases with 
increasing C in the range < ( < 7/27 and expands with C in the range 7/27 < C ^ 1- The 
value Co = 7/27 is due to the half- filled shell (/ = 3) of the up-spin band. Also, we note that 
fs(27, C) > ^f^(oo,C) over the whole range < C ^ 1- Perdew et al. have also studied the 
equilibrium size of spherical clusters of stabilized jellium from a different point of view.El 

Figure |^ shows the ground state E/N and its corresponding for neutral and singly 
ionized jellium clusters as functions of A^, the number of electrons. The energies correspond 
to the equilibrium values for configurations consistent with Hund's first rule. Here, A^ is 
the number of electrons which for neutral clusters equals the number of the ions, and for 
singly ionized clusters it is one less than the number of the ions. The ionization energy of 
an A^-atom jellium cluster is given by 

I{N) = E'°-{N -l,Co)- ^"'"*(iV, Co), (10) 

where ^'^'^"*(A^, Co) is the ground-state energy of the neutral A^-electron jellium cluster which 
occurs at polarization Co consistent with Hund's rule, and E^'^^{N — 1, Co) is the ground state 
energy of the singly ionized cluster with A^ — 1 electrons and A^ positive ions. Obviously, the 
polarizations Co and Co are not the same. In Fig. |^ we have plotted the ionization energies of 
the jellium clusters (A^ < 42) obtained from Eq. (p!0|), and compared them with the results 
obtained from the simple JM calculations (with = 3.99) and experiment on Na clusters. 
The relevance of this comparison with Na results is due to the fact that the equilibrium 
Tg value of the bulk jellium, 4.18, is very close to that of the bulk Na. Except for very 
small clusters, the results based on the equilibrium values show a better agreement with 
experiment than the results of simple JM calculations for Na clusters. 

We saw, in the second step, that if the jellium clusters assume their individual equilibrium 
volumes, then the ionization energies improve. Therefore, we are naturally led to define a 
new set of Ar^ for our SSPJM calculations which are obtained from 

Ar,{N,() = rs{N,()-ao. (11) 

Here, ao = 4.28 for neutral and ao = 4.62 for singly ionized clusters; fs{N, () is the equilib- 
rium Vg value for a cluster with given values of A^ and C- Using Eq. (|ll]), we have performed 
the SSPJM calculations for Na. The results of calculations show that here (in contrast to the 
MSG rule) the energy of a cluster is minimized for a configuration consistent with Hund's 
first rule. This argument is valid for both neutral and singly ionized clusters. This behavior 
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has its roots in the decreasing behavior of fs{N, () over the range < ( < (o for open-shell 
clusters. We have also checked whether or not taking account of the volume change as a 
function of polarization, as given by Eq. (|I]), in the simple JM calculations gives rise to 
the MSG rule. The results show that Hund's first rule remains at work and therefore we 
conclude that, in order to obtain MSG configuration for the ground state, not only Ar^ 
should be an increasing function of C but one should also include the two corrections, due 
to the stabilization, in the simple JM energy. That is, one should use the SSPJM energy 
along with an increasing function for Ar^. This condition will be met if one relaxes the 
spherical constraint on the jellium and considers ellipsoidal shapes for open-shell clusters. 
It is then possible to obtain fs{N,() for a given set of values of N and (, by finding the 
values of the ellipsoid axes that minimize the total energy. We expect that under these 
conditions fs{N,() becomes an increasing function of ( for all clusters. This expectation 
is due to the fact that addition of a nonsymmetric perturbation to a spherical potential in 
a single-particle Hamiltonian lifts the orbital degeneracies and each degenerate level splits 
into (2/ + 1) new levels with different energies. Each of these levels will contain at most 
two electrons with opposite spins. Then, any increment in the polarization is accompanied 
by a transition of a spin-down electron to an unoccupied level and a successive spin-flip. 
This is the only way that one can increase the polarization consistent with Pauli's exclusion 
principle. This process resembles the process of increasing the polarization in a closed-shell 
spherical cluster which results in increasing its equilibrium value. It is a well-known fact 
that the open-shell clusters lose their spherical geometry due to the Jahn-Teller effect,!^ and 
that is why the MSG rule and the odd-even alternations are observed in experimental data 
for alkali metal clusters. In short, using a new set of increasing functions fs{N, () - which is 
obtained using ellipsoidal geometries for simple JM and LSDA- in the SSPJM calculations 
for ellipsoidal geometries of the jellium will lead to the MSG configuration for the ground 
state of the cluster. Work in this direction is in progress. 

Finally, it should be mentioned that in the context of the SSPJM one could fix, at the 
beginning, the pseudopotential core radius for a bulk system of a given species with a given 
fixed polarization; and then use this value in the SSPJM energy functional for a finite cluster 
(transferability condition on the pseudopotential) . El Then, one could obtain the equilibrium 
radius of the jellium by finding that value which minimizes the total energy. We have 
performed such calculations and obtained further agreement to the experimental results, 
which will appear elsewhere. 

IV. CONCLUSION 

In this paper, keeping the spherical geometry, we have considered the finite-size effects 
on the equilibrium values. Our calculations show that for a given A^-electron cluster, the 
quantity fs{N, Q behaves differently for an open-shell and a closed-shell cluster. That is, this 
equilibrium value is an increasing function of ( over the whole range < C < 1 for a closed- 
shell cluster; whereas for an open-shell cluster it is a decreasing function over < ( < (o 
and an increasing function over (o < ( < 1. Here, Co is a polarization corresponding 
to a configuration consistent with Hund's first rule. Our SSPJM calculations based on 
equilibrium values show that, in contrast to the MSG rule, Hund's first rule is at work. 
This behavior is due to the fact that fs{N, Q has a decreasing part for an open-shell cluster. 
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We therefore conclude that to reahze the MSG rule in the SSPJM calculations with fs{N, C), 
and thereby the odd-even alternation, one should lift the spherical constraint on the jellium 
and let it assume ellipsoidal shapes. 
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FIGURES 



FIG. 1. The equilibrium r<j values in atomic units as functions of the polarization Large 
squares, diamonds, crosses, pluses, small squares, triangles correspond to clusters with N=2, 8, 18, 
20, 34, 40 electrons, respectively, for (a) neutral, and (b) singly ionized closed-shell clusters. The 
solid lines and the dashed lines correspond to fitted and bulk values, respectively. 



FIG. 2. The increments of in atomic units as functions of C. The dotted and dashed lines 
correspond to the fitted result [ see Eq. (^) ] for neutral and singly ionized clusters, respectively. 
The solid line corresponds to the bulk jellium. 



FIG. 3. The ionization energies as functions of the number of atoms for Na clusters. The 
diamonds correspond to the dif-SSP JM results obtained using Eq. (|l|) and the triangles correspond 
to the dif-SSP JM results obtained using Eq. (|9|). The small squares are experimental results. 



FIG. 4. The equilibrium values in atomic units and their corresponding total energies per 
electron in electron volts as functions of the number of electrons for neutral clusters with different 
spin configurations. The triangles correspond to the configurations consistent with Hund's first 
rule, and the uppermost symbols for each correspond to the MSG configurations. The symbols 
between the Hund and the MSG configurations correspond to intermediate polarizations. The 
horizontal dashed line below the fg curve corresponds to the value = 4.18 of the bulk jellium. 



FIG. 5. The equilibrium values in atomic units as functions of C- The solid squares and 
diamonds correspond to the closed-shell {N = 34) and open-shell {N = 27) clusters, respectively. 
The dashed line corresponds to the bulk which is obtained using Eq. (17) of Ref. 23. As is seen, 
the equilibrium of the open-shell cluster has a decreasing part before Co = 7/27, whereas for the 
closed-shell cluster it is an increasing function over the whole range < C ^ 1- 



FIG. 6. The equilibrium values in atomic units and their corresponding total energies per 
electron in electron volts for the ground-state configurations (i.e., the configurations consistent 
with Hund's first rule) of A^-electron clusters. The solid squares correspond to neutral and the 
diamonds correspond to singly ionized A'^-electron clusters. As is expected, for small clusters the 
differences are high and for large clusters these differences approach zero. 



FIG. 7. The ionization energies in electron volts. The diamonds correspond to the simple JM 
results for Na using the bulk value = 3.99. The small solid squares correspond to the ionization 
energies of jellium clusters in their equilibrium states obtained using Eq. (p^. It shows a good 
agreement with the experimental values of Na, which are shown as the large squares. 
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